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Abstract 

We present model independent measurements of the helicity basis form factors in the decay 

D'^ K^TT^e'^Uc obtained from about 2 800 decays reconstructed from a 281 pb~"^ data sample 
collected at the V'(3770) center-of-mass energy with the CLEO-c detector. Wc confirm the existence 
of a previously observed spin-zero K~7r~^ component interfering with the K amplitude. We see 
no evidence for additional d- or /-wave contributions. 
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I. INTRODUCTION 



Exclusive semileptonic decays are excellent probes of charm decay dynamics since strong 
interaction effects only enter through the current coupling the parent D meson with the final 
state hadronic system This current is generally expressed in terms of form factors that 
depend on g^, the squared invariant mass of the virtual materializing as a positron- 
neutrino pair. The specific dependence is an ansatz in most of the models, although 
the ultimate goal is to derive both normalization and dependence from first principles. 
Various theoretical approaches have been pursued to study semileptonic form factors, such 
as quark models ^5], QCD sum rules j^, lattice QCD and the parameterization of the 
form factor dependence based on effective poles jsj constrained by heavy quark and soft 
collinear effective theories. Measurements of both the branching fraction |8|| and the q^ 
dependence of the form factors have been made for the — > K*^l'^i'i decays using 
specific parameterizations Recently, Ref. reported a first observation of an additional 
component besides the dominant K £"'"z/£ in the K^n'^i'^ui decay process. 

In the additional component, K~ and vr^ are in a relative s-wave, which has revealed an 
interesting connection between semileptonic decays and light quark physics. According to 
Ref. 0, 2.4 ± 0.7% of the decays in the mass range 0.8 GeV/c^ < mK-K < 1-0 GeV/c^ are 
due to this s-wave component, where rnxn is the K~7c^ mass. 

Using a technique developed by FOCUS [l^, we present non-parametric measurements 
of the dependence of the helicity basis form factors that describe an amplitude for the 
K~n~^ system to be in any one of its possible angular momentum states. This is done 
by projecting out the helicity form factors directly from data without the use of fitting 
functions. Our results will allow theorists to directly compare their models with the data 
free from parameterization. We also confirm the existence of the s-wave component, study its 
form factor for the first time, and limit the strength of possible d- and /-wave contributions. 

There are several good reasons to perform this analysis using the CLEO-c data: The 

resolution in CLEO-c is an order of magnitude better than in FOCUS. In addition, the 
K e^Vf, process is simpler than K //"''z/^, since there is a negligible probability 

for the much less massive to be left-handed which eliminates one of the form factors. The 
additional form factor describing the left-handed coupling has an angular distribution so 
similar to Ho in Eq. (P) that it degrades the measurement of Hq by a factor of three. 

The amplitude Ai for the semileptonic decay K^n^e^Ve is described by five 

kinematic quantities: g^, mK-n, the angle between the vr and the D direction in the K~7i^ 
rest frame {Ov), the angle between the and the D direction in the e+z/g rest frame (6'^), and 
the acoplanarity angle between the two decay planes (x)- Following Ref. Q, we can express 
the matrix element for the decay — > K^n^e'^Ug in the vicinity of the K mass in terms 
of the three hehcity amplitudes (if+(g^),if_(g^), and Ho{q^)) describing the pseudoscalar to 
vector hadronic current and one form factor {hQ{q^)) describing a broad s-wave resonance. 
The dominant terms in the differential width integrated over the angle x can be 

expressed as: 

j \M\''dxo^Gl\V,s\\q^-ml) 

X [ ((1 -|- cos Oi) sin 

+ ((1 - cos^£) sin6vf\H^{q^)\^\BW\^ (1) 
+ (2sin^^cosey)^|ifo(g^)n5W^P 
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+8 sin^ 9i cos 9vHo{q'^)ho{q^) Re{Ae-'^BW} 

The second to last term in Eq. represents the interference between the s-wave K~tt^ and 
the K amphtude, where A and S represent the amphtude and the phase of the s-wave, 
respectively. Since the s-wave amplitude could be observed only through interference with 
the K , FOCUS was sensitive to the s-wave amplitude only in the vicinity of the K pole. 
They modeled the s-wave contribution as a constant complex amplitude, even though a 
variation of the amplitude with respect to rriKn is possible. In this paper, we assume the A 
and 6 values obtained in Ref. 0] in the study of the hQ^q^) HQ^q^) interference term. The 



K amplitude in Eq. ([T]) is represented as a Breit-Wigner of the form: 

BW 
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where P* is the kaon momentum in the K~t[^ rest frame and Pg* is the kaon momentum in 
this frame at the resonant mass mo. The x integration significantly simplifies the intensity 
by eliminating all interference terms between different helicity states of the virtual with 
relatively little loss in form factor information. We have dropped the term which is second 
order in the small, s-wave amplitude A. 

The three helicity basis form factors for the K e+z/g component are generally 

written jlj] as linear combinations of vector {V^q"^)) and axial- vector {Ai^2{q^)) form factors 
according to Eq. Q: 



H^q') = (Mb + ruKMiiq') T 2 ^^^^^ Viq') and (2) 



Ho{q 



2\ 



2mx7rv^ 



{Ml - mL - q'){Mo + rriKMAq') - A,{q') 



where M£, is the mass of the and K is the momentum of the i^~7r+ system in the rest 
frame of the . 

In this paper, we use a projective weighting technique ^3 to disentangle and directly 
measure the q^ dependence of these helicity basis form factors free from parameterization. 
We provide information on the form factor products HQ^q^), and ho{q'^) Holq'^) in 

bins of by projecting out the associated angular factors given by Eq. ([T)). 



II. EXPERIMENTAL AND ANALYSIS DETAILS 

The CLEO-c detector 2| consists of a six-layer, low-mass, stereo-wire drift chamber, 
a 47- layer central drift chamber, a ring-imaging Cerenkov detector (RICH), and a cesium 
iodide electromagnetic calorimeter inside a superconducting solenoidal magnet providing a 
1.0 T magnetic field. The tracking chambers and the electromagnetic calorimeter cover 93% 
of the full solid angle. The solid angle coverage for the RICH detector is 80% of Air. Iden- 
tification of the charged pions and kaons is based on measurements of specific ionization 
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FIG. 1: The niK-K distribution for events satisfying our nominal D'^ — > K~TT^e^i'e selection 
requirements. Over the fuh displayed mass range, there are 2 838 events satisfying our nominal 
K^TT^e^Ve selection. For this analysis, we use a restricted mass range 0.8 — 1.0 GeV/c^ 
(shown by the vertical dashed lines). In this restricted region, there are 2 472 events. 




(dE/dx) in the main drift chamber and RICH information. In positron identification, in ad- 
dition to dE/dx and RICH information, the ratio of energy deposited in the electromagnetic 
calorimeter to the measured track momentum (E/p) is used. 

In this paper, we use 281 pb^^ of data taken at the '?/'(3770) center-of-mass energy with 
the CLEO-c detector at the Cornell Electron Storage Ring (CESR), which corresponds to 
a sample of 0.8 million D^D~ pair events 113 . Monte Carlo (MC]events are generated by 
EVTGEN [l3| and the detector is simulated using a GEANT-based |lfi| program. Simulation 
of final state radiation is handled by photos|]J|. Throughout this paper charge-conjugate 
modes are implied. 

We select events where a —>■ K~Tt^e^Vf. candidate is produced against a fully recon- 
structed tagging D~ . The tagging decays into one of the following six decay modes: 
D- K^TT-, D- K+n-n-, K^ti-tt^, ^ K+'r-ti-'k^, ^ K^g^^- n- n+ , and 

K^K^n^ . Multiple D~ candidates per event are allowed. More details on selecting 
the t agg ing D~ candidates as well as identifying 7r° and Kg candidates are described in 
Ref. [if. 

The following selection cuts represent our nominal selection criteria. The semileptonic 
reconstruction starts by requiring three well-measured tracks not coming from the tagging 
decay in the event. Charged kaons and pions are required to have momenta of at least 
50 MeV/c and are identified using dE/dx and RICH information. Positron candidates are 
required to have momenta of at least 200 MeV/c, satisfy | cos6'| < 0.9, where 6 is the angle 
between the positron and the beam line, and pass a requirement on a likelihood variable 
that combines E/p, dE/dx, and RICH information. We obtain 2 838 K^Ti^e^Ug 
events. Finally, we require 0.8 GeV/c^ < < 1-0 GeV/c^ and select 2 472 events. The 
rnxn distribution for these K^n^e^Uf, candidates is shown in Fig. ^ In our analysis, 

the only particle not reconstructed is the neutrino of the K'lr^e'^Vf. decay. The 

neutrino four-momentum can be reconstructed from the energy-momentum balance in the 
beam-beam center-of-mass frame. We assign the neutrino the missing energy (ii^miss) which 
is the difference between the beam energy (-Ebeam) and the sum of the energies of the charged 
semileptonic daughters. The neutrino momentum is obtained from direction of the missing 
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FIG. 2: The observed energy imbalance U distribution for K^ii^e'^Ve events satisfying our 

nominal selection. We show the boundaries for the additional \U\ < 0.06 GeV requirement used 
for the tighter selection. The inset is the same distribution on a finer scale where the background 
contribution estimated from a Monte Carlo sample is more apparent as a hashed histogram. 

momentum (p^iss = Pmiss/|Pmiss|) defined as: 

Pmiss = -\[^L~^rr^pD- - PK- - P7T+ - Pe+- (3) 

In Eq. (jni), Pd- is the measured direction of the tagging D~ , pis the measured momentum 
of the indicated particle, and m^) is the known mass of the D~ fv?\. The magnitude of the 
neutrino momentum is set to -Emiss- 

We will compare the nominal selection to a tighter selection where we add two further 
restrictions: 0.846 GeV/c^ < niK^ < 0.946 GeV/c^ and \U\ < 0.06 GeV, where U = 
Emiss — c|pmiss|- Figure 12 shows the U distribution for the nominally selected sample with 
the background contribution. The background shape in Fig. |21 is obtained using a charm 
MC set which consists of generic DD events. We estimate a background level of 3.6% for the 
nominal selection and 0.4% for the tighter selection. We have checked the MC background 
simulation with a sample of wrong sign events where we reconstructed the semileptonic 
candidates with electrons instead of positrons and obtained consistent results between the 
data and the MC samples. 

The resolution predicted by our MC is roughly Gaussian with an r.m.s. width of 0.02 
GeV^/c'^, which is shown in Fig. 01 This is negligible on the scale in which we will bin our 
data. 

Even though we only measure the shape and not absolute rates, as a cross-check we have 
measured the K e^Ve branching fraction assuming an efficiency derived from the 

generic MC. We found statistical consistency with the result in Ref. [18ii |. 

III. PROJECTION WEIGHTING TECHNIQUE 

We extract the helicity basis form factors using the projective weighting technique more 
fully described in Ref. [lO|. For a given bin, a weight is assigned to the event depending 
on its Oy and 6^ decay angles. We use 25 joint Acos6'y x ^cosOn angular bins: 5 evenly 
spaced bins in cos6'v times 5 bins in cosO^. 
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FIG. 3: The resolution predicted by the generic DD MC. 

For each qf bin, the angular distribution can be written as a vector Ni whose components 
give the number of events reconstructed in each of the 25 angular bins. According to Eq. (Q), 
Ni can be written as a linear combination of vectors with coefficients faiQi), 

M = /+ (gf) m+ + /_ (gf) m_ + /o (q^) mo + // (g- ) mj . (4) 

The m_|_, m_, mo, and rnj vectors are computed using a full detector simulation weighted by 
the corresponding helicity term in Eq. (^. For example, m+ is computed using a simulation 
generated with an arbitrary function for (such as = 1) and zero for the 

remaining three form factors. The fa{(li) functions are proportional to the true H^{qf) along 
with multiplicative factors such as G"^ \ Vcsf {q^ """^f) and additional corrections accounting 
for experimental effects such as acceptance. 

Ref. shows how Eq. Q can be solved for the four form factor products if^(g^), 
if^(g^), Hl{q^), and hQ{q^)HQ{q^) by making four weighted q^ histograms. The weights are 
directly constructed from the four fha vectors. We have tested the method using thirty 
independent Monte Carlo samples. Figure |3] demonstrates that the projective weighting 
technique returns realistic errors with no significant bias. 

IV. RESULTS 

Figure El shows the four form factor products multiplied by q^ obtained for data using 
the nominal selection criteria. In this Figure, the background estimated from the generic 
charm MC simulation is subtracted. For the nominal selection we estimate that 3.6% of 
the events are background. The form factors are normalized by scaling the four weighted 
histograms (if^(g^), Hl{q^), and hQ{q^)HQ{q^)) by a common factor so that q^Hl{q^) = 1 as 
q^ 0} Because of our excellent q^ resolution, there is negligible correlation among the six 
q^ bins for a given form factor product, but the relative correlations between different form 



The details of the normalization procedure are as follows: The solid curve in Fig. [3 (c) represents the 
parameterized model for the form factor q^HQ^q^) presented in Ref. 01 ■ This curve is normalized to 1 at 
q^ ^ Q. Then our data points in Fig.[Sl(c) are scaled to the curve by a factor obtained from a fit. The 
data points in the other Fig. histograms are all scaled by the same factor. 
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FIG. 4: Monte Carlo study of the projective weighting technique errors and biases. The recon- 
structed form factors are the points with error bars. The curves give the hehcity form factor 
products assumed in the simulations. Thirty independent simulations were used for each bin. 
The averages of the simulated results are given by the horizontal lines. Their actual ordinate is at 
the left of the averaging line. The yield in these MC samples are six times the yield in the data. 
The plots are: (a) Hl{q^), (b) Hl{q^), (c) H^q^), and (d) ho{q^) Ho{q^) . 



factor products in the same bin can be as high as 36%. The /io(Q'^) HQ{q^) form factor 
in Fig. El (d) is measured through the interference term of Eq. (0), which is proportional to 
Re{Ae-'^BW}. Avera ging over the full trk-k mass range, Re{^e BW} is proportional to 
—Asm 5. The interference term is then measured by a fourth projective weight and divided 
by the FOCUS —A sin 5 value to obtain ho{q^) Ho{q'^). The measurement on ho{q'^) Ho{q'^) 
includes a systematic uncertainty due to uncertainties in A sin 6, which we discuss below. 
Figure El shows the helicity amplitudes without weighting, thus emphasizing the low 
region. It includes results obtained with the nominal selection criteria as well as the more 
restrictive selection criteria. The background fraction in the tightly selected sample is only 
0.4% of the event sample and its background subtraction is negligible. The consistency of 
the two results shows that the systematic error due to background subtraction is small. 

We have considered several sources of systematic uncertainties. Even though we have very 
little background in our data set, the background subtraction still constitutes the primary 
source of systematic error. For this uncertainty, we assign a conservative value by reducing 
the level of the background being subtracted by a factor of two and comparing these form 
factor products with the results based on the full background subtraction. For if^(g^), the 
background systematic uncertainty is estimated to be 16.5% for the first bin and less than 
4.8% for higher bins. For i7i(g^), the background systematic uncertainty is estimated to be 
41.3% for the first g^ bin and less than 8.1% for higher bins. For ifQ(g^), the background 
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FIG. 5: Measured form factor products for six bins. The data are points with error bars 
representing only the statistical uncertainties. The solid curves are obtained using the form 
factor representation and the parameters reported in Ref. 01, namely: V{Q)/Ai{Q) = 1.505, 
^2(0)7^1(0) = 0.875, and s-wave parameters A = 0.33 and 5 = 39°. The histogram plots are: 
(a) q^Hl(q^\ (b) q^Hl{q^), (c) q^Hl(q^\ and (d) q^K{q^)Ho[q^). 
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FIG. 6: Measured form factor products for six q^ bins. We show the results obtained with the 
nominal selection criteria (1st point in each bin) and the tighter selection criteria (2nd point in 
each bin). The error bars represent statistical uncertainties only. The solid curves are obtained 
using the form factor representation and the parameters reported in Ref. Q. The histogram plots 
are: (a) (g^), (b) Hliq'), (c) HHq'), and (d) ho{q')Hoiq^). 



systematic uncertainty is estimated to be less than 2.1%; and for hoHo, the background 
systematic uncertainty is estimated as less than 18.0%. The systematic uncertainty due 
to track reconstruction and particle identification is rather low since we are reporting on 
form factor shapes rather than absolute normalization. The systematic uncertainty from 
these sources is estimated as less than 1.9% for all of the form factor products. Lastly, 
we assess an additional scale error of 13.4% for the hnHo form factor product due to the 
uncertainties in the A and 6 values reported in Ref. |7|. Our total systematic uncertainty 
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is the quadrature sum of these systematic contributions and found to be small compared to 
the statistical uncertainty. In all cases except hoHo, the systematic uncertainty is dominated 
by the background systematic uncertainty. Table H] gives the range for each bin, along 
with the form factor products, their statistical uncertainty (first error) and the estimated 
systematic uncertainty (second error). 



TABLE I: Summary of form factor product results for six q'^ bins. Each form factor product is 
averaged over the indicated range. The first and second errors are statistical and systematical 
uncertainties, respectively. The numbers are normalized using the condition: q'^H^lq'^) = 1 as 
g2 ^ 0. 

range (GeVVc^) Hi hpHo 

0.000 - 0.167 0.16 ± 0.12 ± 0.03 1.19 ± 0.87 ± 0.49 13.17 ± 0.91 ± 0.13 13.60 ± 4.72 ± 1.86 
0.167 - 0.334 0.32 ± 0.10 ± 0.02 1.99 ± 0.34 ± 0.16 4.30 ± 0.34 ± 0.04 6.73 ± 1.84 ± 0.92 
0.334 - 0.501 0.52 ± 0.11 ± 0.01 1.88 ± 0.23 ± 0.08 2.59 ± 0.23 ± 0.04 3.34 ± 1.25 ± 0.46 
0.501 - 0.668 0.61 ± 0.10 ± 0.01 1.93 ± 0.17 ± 0.08 1.69 ± 0.18 ± 0.01 -0.36 ± 0.94 ± 0.08 
0.668 - 0.835 0.55 ± 0.10 ± 0.01 1.53 ± 0.15 ± 0.03 1.43 ± 0.18 ± 0.03 0.35 ± 0.87 ± 0.06 
0.835 - 1.000 0.65 ± 0.18 ± 0.01 1.46 ± 0.23 ± 0.03 1.47 ± 0.25 ± 0.02 -1.78 ± 0.84 ± 0.24 



Figures El (d) and El (d) show that the ho{q'^)HQ{q'^) form factor product is in rough 

agreement with the model from Ref. 0]. We are also consistent with the FOCUS s-wave 

phase (5 = 39 ± 5°. Our consistency check is the comparison of the s-wave interference term 

for events with mK-K below and above the K pole as shown in Fig. [71 The interference 

between the s-wave and the Breit-Wigner is proportional to cos{S — (3±) where /?_ and /3+ 

— *o 

are the average Breit-Wigner phases below and above the K pole, respectively. In the 
mass region where m{K ) < niKn < 1-0 GeV, /?+ is —59°, which is nearly orthogonal to 
the FOCUS s-wave phase. If our phase is consistent with the FOCUS phase, the effective 
ho{q'^) HQ{q'^)Re{Ae~'^^ (BW)} should disappear above the K*^ pole as it does in Fig. [71 (b). 
The sum of the first three bins of the histogram of Fig. |71 (a) differs from zero by 5.9 
standard deviations. This confirms the earlier observation of an s-wave interference with 

-ho(q^)Ho(q^)Re{/\e-(BW)} 

I " " I " " j 
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FIG. 7: The s-wave interference term for the events (a) below the K pole and (b) above the 
pole. The interference term depends on the s-wave phase relative to the average phase of the 
Breit-Wigner in the corresponding plot. All of the cosOy interference observed by FOCUS ^ was 
also below the K pole. 

the FOCUS phase, but does not provide an independent measurement of that phase. 
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It is of interest to search for the possible existence of additional non-resonant ampli- 
tudes of higher angular momentum. It is fairly simple to extend Eq. (Q) to account for 

potential ci-wave or /-wave interference with the K Breit-Wigner. We search specif- 
ically for a possible zero helicity d-wave or /-wave piece that interferes with the zero 



-*o 



helicity K"" contribution. One expects that potential hQ^{q^) and hQ \q^) form factors 
would peak as 1/ \/q^ near ^ as is the case for the other zero helicity contributions 
Ho{q'^) and hQ{q^). If so, the zero helicity contributions should be much larger than po- 
tential d- or /-wave ±1 helicity contributions. The d- and /-wave interference terms are 
the same as the s-wave interference term apart from additional Wigner (i-matrices. These 
describe anisotropy in the d- and /-wave decays. The additional factors are (^oo(^v') for 
the d-wave and d^Q{6v) for the /-wave. Hence, the (i-wave weights are based on an fifth 

term of the form 4 sin^ 0^ cos Oy (3 cos^ Oy - 1) HQ{q^) hf{q^) Re{Ad e-'^-^BW} in Eq. jH). 
The /-wave weights are based on an fifth term of the form 4 sin^ Oi cos 9v (5 cos^ 9v — 
3 cos 6v) HQ{q'^) h'Q\q'^)'Re{Af 6'"^^^ BW}. Averaging over the Breit-Wigner, the interfer- 
ence should be proportional to Adj sm6dj h^o'^\q^) Ho{q^)- 
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FIG. 8: Search for (a) d-wave and (b) /-wave interference effects, Adj sin 6dj h'^'^\q'^) HQ{q'^), as 
described in the text. 



Figure IHl shows the Adj sin 6dj hlf''^\q'^) HQ{q'^) obtained with this technique. It is ap- 
parent from Fig. |H1 that we have no compelling evidence for either a ci-wave or an /- 
wave component. Unfortunately, we cannot set a particularly stringent limit. Under 
the assumption h^Q'^\q^) = Holq"^), used in Ref. Q, we perform a fit of Fig. |H1 to 
the form Adj sin 5dj H^^q^). The results of these fit were Ad sin5d = —0.07 ± 0.20 and 
Af sin 5/ = 0.17 ± 0.18. For comparison, the value for A sin 5 for the s-wave contribution 
according to Ref. is 0.21 ± 0.028. 



V. SUMMARY 

We present a non-parametric analysis of the helicity basis form factors that control the 
kinematics of the K'n^e^Ug decays. We use a projective weighting technique that 

allows one to determine the helicity form factor products by weighted histograms rather 
than likelihood based fits. The resulting helicity form factors have very little background 
and the background subtraction procedure is well understood. We confirm the existence 
of an s-wave K~tt^ component to the K~Ti^e^Ve decay and have studied the q^ 
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dependence of its form factor. Finally, we have searched for possible d-wave or /-wave non- 
resonant interference contributions to K^n^e^Ve- Although we have no statistically 
significant evidence for f-wsive or (i-wave interference, we are unable to limit these terms to 
appreciably less than the established s-wave interference. 
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Foundation. This work was supported by the National Science Foundation, the U.S. De- 
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